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1 Introduction 

One of the central mysteries of quantum theory is that it seems to be fundamen- 
tally nonlocal - and yet the nonlocality cannot be used for practical signalling 
at a distance. As argued elsewhere (Valentini 1991a, b; 1992; 1996), the consis- 
tency of modern physics seems to depend on a 'conspiracy', in which nonlocality 
is hidden by quantum noise. It is as if there is an underlying nonlocality which 
we are unable to control because of the statistical character of quantum events. 

A natural explanation for this peculiar state of affairs is the hypothesis that 
quantum probability distributions are not fundamental, but merely represent a 
special state of equilibrium in which nonlocality happens to be masked by sta- 
tistical noise (Valentini 1991a, b; 1992; 1996). In quantum theory a system with 
wavefunction tp has an associated probability distribution given by the Born rule 
p — \ip\ 2 . This is usually regarded as a fundamental law. If instead we regard 
p — |'0| 2 as analogous to, say, the Maxwell distribution of molecular speeds for a 
classical gas in thermal equilibrium, the above 'conspiracy' takes on a different 
light: it is seen to be an accidental, contingent feature of the 'quantum equilib- 
rium' distribution p — \ip\ 2 . In a universe that is in (global) thermal equilibrium, 
it is impossible to convert heat into work; this is the classical thermodynamic 
heat death. In a universe that is everywhere in quantum equilibrium, it is im- 
possible to use nonlocality for signalling; this is the 'subquantum heat death' 
which, we claim, has actually occurred in our universe. 

The pilot- wave theory of de Broglie and Bohm (de Broglie 1928; Bohm 
1952a, b) offers a concrete model of this scenario.^] A system with wavefunction 
tjj is assumed to have a definite configuration x(t) at all times whose velocity 
is determined by the de Broglie guidance equation n = dS/dx, where 7r is the 
canonical momentum and S is the phase of ip (given locally by S — filming). 
The wavefunction ip is interpreted as an objective 'guiding field' in configuration 
space, and satisfies the usual Schrodinger equation^] 

0a Present address: Theoretical Physics Group, Blackett Laboratory, Imperial College, 
London SW7 2BZ, England, and Center for Gravitational Physics and Geometry, Department 
of Physics, The Pennsylvania State University, University Park, PA 16802, USA. Permanent 
address: Augustus College, 14 Augustus Road, London SW19 6LN, England. 

iNote that in 1927 at the Fifth Solvay Congress de Broglie proposed the full pilot- wave 
dynamics in configuration space for a many-body system, not just the one-body theory; and 
unlike Bohm's reformulation of 1952, de Broglie's new approach to dynamics in the 1920s was 
always based on velocities (Valentini 2001a). The proceedings of the Fifth Solvay Congress 
are being translated into English (Bacciagaluppi and Valentini 2001). 

2 Pilot- wave dynamics may be applied to fields as well as particles — even to fermion fields 
(Valentini 1992; 1996; 2001a) and non-Abelian gauge theories (Valentini 2001a). Note that 
the natural spacetime structure associated with this velocity-based dynamics is Aristotelian 
spacetime E X _E 3 , which has an inbuilt natural state of rest (Valentini 1997). 
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At the fundamental hidden-variable level, pilot-wave theory is nonlocal. For 
example, for two entangled particles A and B the wavefunction ip(xA , xb , t) has a 
non-separable phase S(xa,xb, t) and the velocity dxA/dt = V AS(xA,XB,t) /m 
of particle A depends instantaneously on the position xb of particle B.^| How- 
ever at the quantum level, where one considers an ensemble with distribution 
p(xA,XB,t) = \iP(xa, xb, t)\ 2 , operations at B have no statistical effect at A: as 
is well known, quantum entanglement cannot be used for signalling at a distance. 

It is worth emphasising that this 'washing out' of nonlocality by statistical 
noise is peculiar to the distribution p = \tp\ 2 . If one considers an ensemble of 
entangled particles at t = with distribution po [xa , Xb ) ^ \tI>o(xa, xb)\ 2 , it may 
be shown by explicit calculation that changing the Hamiltonian at B induces an 
instantaneous change in the marginal distribution pa(xa, t) = JdxB p(xa, xs,t) 
at A. For a specific example it was found that a sudden change Hb — > H' B at 
B - say a change in potential - leads after a short time e to a change ApA = 
Pa(xa,c) — Pa{xa,0) at A given by (Valentini 1991b) 

A e 2 9 ( , . f ' . po{x A ,x B ) - \ip {x A ,x B ) | 2 \ , . 

Ap A = - — - — [a(x A ) / dx B b{x B ) — — (1) 

AmdxA V J \ipo{XA,x B )r J 

(Here cl{xa) depends on ipo; b{xs) also depends on H' B and is zero if H' B = Hb-) 
For nonequilibrium ensembles po ^ IV'ol 2 ; there are genuine instantaneous 
signals at the statistical level. This has recently been shown to be a general 
feature, independent of pilot-wave theory. Any deterministic hidden-variables 
theory that reproduces quantum theory leads, for hypothetical nonequilibrium 
ensembles, to instantaneous signals at the statistical level (Valentini 2001b). 

Note that our inability to see the trajectories - or hidden variables - is 
also a contingent feature of equilibrium: the uncertainty principle holds if and 
only if p — \ip\ 2 (Valentini 1991b)J^| Heuristically, it is natural to compare our 
limitations with those of a Maxwell demon in thermal equilibrium with a gas, 
whose attempts to sort fast and slow molecules fail. The common objection to 
hidden variables - that their detailed behaviour can never be observed, making 
their existence doubtful - is seen to be misguided: for the theory cannot be 
blamed if we happen to live in a state of statistical equilibrium that masks the 
underlying details. There is no reason why nonequilibrium could not exist in 
the remote past or in distant regions of the universe (Valentini 1992), in which 
case the details of the 'hidden-variable level' would not be hidden at all. 

Our central theme, then, is the subquantum 'heat death of the universe' 
at the hidden-variable level. How does one account for the quantum noise - 
encapsulated by the Born rule p = \tp\ - that pervades our observed universe 
at the present time? 

3 It is being assumed here that the position of particle B could in some sense be changed 
while keeping the wavefunction fixed. 

4 Note, however, that the approximately classical functioning of the experimenter also plays 
a role here. A subquantum demon - an automaton that functions at the hidden-variable level 
— would see subquantum trajectories even for an equilibrium ensemble (Valentini 2001a). 
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2 'Empirical' Approach to Statistical Mechanics 



It will be shown how one may set up an analogue of classical statistical mechan- 
ics, based on pilot-wave dynamics in configuration space (Valentini 1991a; 1992; 
2001a). We shall focus in particular on an analogue of the classical coarse- 
graining approach. But first, we shall address some foundational issues that 
were not considered in the author's original papers. 

In what sense might one 'explain' the probability distribution p = \tp\ 2 , in 
the context of a deterministic (and time-reversal invariant) pilot-wave theory? 

Say the Moon is now at position P (in phase space) at time t. This might be 
'explained' - using Newton's laws - by the fact that the Moon was (earlier) at 
position Q at time to; and if to is in the remote past, one would in practice deduce 
that the Moon must have been at Q at time to, from the observed position P 
today. That one has had to deduce the past (Q) from the present (P) would 
not change our physical intuition that the Moon may reasonably be said to be 
at P now 'because' it was at Q at to- 

This example may seem remote from statistical mechanics. However, many 
workers would object - in the author's opinion wrongly - that Q at to is a mere 
deduction, and cannot be regarded as a satisfactory explanation for P today. 
The argument would be that one would have a satisfactory explanation only if 
it could be shown that all - or at least 'most' - possible positions of the Moon 
at t evolve into the Moon being at P today. Only then would one have an 
'explanation' for the present, as opposed to a mere deduction about the past. 

If this seems an unfair portrayal of what many workers in the foundations of 
statistical mechanics would claim, then consider what is attempted in that field. 
The 'Holy Grail' has always been to explain relaxation to thermal equilibrium 
by showing that all or 'most' initial states do so. And the great stumbling block 
has always been that, by time-reversal invariance, for every initial state that 
evolves towards equilibrium there is another that evolves away from it. One 
must then place restrictions on initial states, such as an absence of correlations, 
or a lack of fine-grained microstructure, and so on, restrictions that are violated 
by the unwanted 'time- reversed', entropy-decreasing states. One then tries to 
argue that the restricted set of initial states is 'natural', or that 'most' initial 
states - with respect to some suitable measure - satisfy the required conditions. 
As is well known, this program is as controversial today as it was a century ago. 

To avoid such controversy one might adopt a more modest - and more real- 
istic - 'empirical' approach (Valentini 1996; 2001a). Going back to the example 
of the Moon, one would try to deduce - or perhaps guess - where the Moon 
must have been before to explain its position today. In the same sense, for a box 
of gas evolving to thermal equilibrium, one should try to deduce or guess what 
the initial (micro-)state must have been like to yield the observed behaviour - 
without trying (in vain) to show that all or 'most' initial states would do so. Of 
course, in the case of a gas there will be a whole class of microstates yielding 
the required behaviour; one will not attempt to deduce the exact initial state 
uniquely (unlike in the case of the Moon). Further, because there are so many 
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variables one resorts to statistical methods. Thus, one tries to construct a class 
of initial conditions that yields the observed behaviour, and one tries to under- 
stand the evolution of that class towards a unique equilibrium state on the basis 
of a general mechanism (without having to solve the exact equations of motion). 

In the case at hand, what is to be explained is the observation of equilibrium 
today to within a certain experimental accuracy, p = |V>| 2 ±e- For example, for a 
large number of independent Hydrogen atoms each with ground-state wavefunc- 
tion V'ioOj one might measure the distribution of electron positions (with respect 
to the nuclei, in a nonrelativistic model) as accurately as current technology al- 
lows, and obtain a result p = IV'iool 2 ± e. The initial or earlier conditions of the 
universe must have been such as to reproduce this result today. 

One possible initial condition is just equilibrium itself (which is preserved by 
the equations of motion). That is, the universe could have started in equilibrium, 
leading to p = \xp\ 2 exactly today. But initial equilibrium is only one possibility 
among (uncountably) many: even without any proofs or assumptions about 
relaxation p — > |^| 2 , it is clear on grounds of continuity alone that if initial 
equilibrium evolves to equilibrium today then an infinite class of initial states 
sufficiently close to equilibrium must evolve to p = \ip\ 2 ± e today; and given 
the violence of the early universe, one expects there to be a large class of initial 
states far from equilibrium that evolve to p = |?A| 2 ± e today as well. 

What sort of past conditions could have evolved to p = \tp\ 2 ± e today? 
Below, we shall hypothesise a class of possible initial (nonequilibrium) states 
and give a general mechanism explaining their approach to equilibrium, as is 
done in classical statistical mechanics. 

It is not true, of course, that all initial states evolve towards equilibrium; 
time-reversal invariance forbids. Nor does it need to be true: 'bad' initial con- 
ditions are ruled out on empirical, observational grounds. 

3 Subquantum if-Theorem 

There are many approaches to classical statistical mechanics. Here we focus 
on a pilot-wave analogue of the classical coarse-graining H -theorem (Valentini 
1991a). This is not because we believe coarse-graining to be the best approach; 
it has advantages and disadvantages, like the others. 

For a classical isolated system, both the probability density p and the vol- 
ume element dil (on phase space) arc preserved along trajectories. The clas- 
sical //-function H c i ass = J dQ plnp is constant in time. If we replace the 
fine-grained p by the coarse-grained p and assume that po — po at t = 0, then 
H c iass{t) < H c iass{0) - the classical coarse-graining fl-theorem of the Ehren- 
fests. The decrease of Hciass corresponds to the formation of structure in p and 
the consequent approach of p to uniformity. (See for example Tolman (1938).) 

Consider now an ensemble of complicated many-body systems, each with 
wavefunction 'J, the configurations X distributed with a probability density P. 
The gradient V S of the phase of \& determines a velocity field X in configuration 
space. (For a low-energy system of N particles of mass m, we have 3-vector 
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velocities X{ = ViS/m where i = 1, 2, N.) The continuity equations 

| + V'(XP)=0 (2) 



(which holds by definition of P) and 

d\V\ 2 



dt 



V-(X|*| 2 ) = (3) 



(which follows from the Schrodinger equation) imply that the ratio f = P/\^f\ 2 
is preserved along trajectories: df /dt — 0, where d/dt = d/dt + X ■ V. Further, 
like dfl classically, the element |^| dX is preserved along trajectories too (where 
dX is the volume element in configuration space). This suggests replacing p — > / 
and d£l — > |\I/| 2 (iV in H class , yielding the subquantum H- function 

H = j dX Pln(P/|*| 2 ) (4) 

This is in fact the relative negentropy of P with respect to l^l 2 . 

The above continuity equations now imply that H = H(t) is constant: 
dH/dt = 0, as in the classical case. But if one divides configuration space 
into cells of volume SV and averages P and l^l 2 over the cells, yielding coarse- 
grained values P and |W| 2 , one may define a coarse-grained //-function 

H = j dX P\n(P/\^) (5) 

If we assume the initial state has 'no fine-grained microstructure' at t = 0, 



Po = Pa, |*o| 2 = |*or (6) 

then it may be shown that 

Ho - H(t) = jdX\V\ 2 (/ ln(///) + / - /) > (7) 

(where / = P/WJ 2 ), so that 

H(t) < H (8) 

for all t (Valentini 1991a). 

It is instructive to examine the time derivatives of H(t) at t = 0. Use of the 
continuity equations shows that (dH/dt) Q = and (Valentini 1992) 

(^f ) = " / ^ - * ° (9) 

The quantity in brackets is just the (non-negative) variance var<5y(X • V/o) of 
(Xq- V/o) over a coarse-graining cell SV. If we suppose that var^y (-^o • V/o) ^ 
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then we have the strict inequality H(t) < Hq and H(t) must decrease over at 
least a finite time interval (0,T) immediately after t — (Valentini 2001a). 

Because xhi(x/y) + y — x > for all real x, y, the coarse-grained //-function 
is bounded below by zero, H > 0; further, H — if and only if P = \fy\ 2 
everywhere (Valentini 1991a). A decrease of H towards its minimum value then 
corresponds to an approach of P towards |\&| 2 . 

The decrease of H corresponds to a 'stirring' of the two 'fluids' P and l^l 2 
by the same velocity field X, making P and |^| 2 less distinguishable on a coarse- 
grained level. (This is similar to the classical Gibbs stirring of two liquids.) As 
in the corresponding classical case, the //-theorem gives us an insight into the 
mechanism whereby equilibrium is approached (for the particular class of initial 
conditions specified by (6)). Whether or not equilibrium is actually reached 
will depend on the system. It must be assumed that the initial steps towards 
equilibrium described by the //-theorem are actually completed in Nature, for 
appropriately complex systems, so that H(t) — > and P — ► |^| 2 . 

Given (coarse-grained) equilibrium P = \^\ 2 for our ensemble of many- 
body systems, it is straightforward to show that if a single particle is extracted 
from each system and prepared with wavefunction ip, the resulting ensemble of 
particle positions has a (coarse-grained) distribution p = \ip\ 2 (Valentini 1991a). 
It is therefore to be envisaged that the equilibrium distribution seen today 
for ensembles of single particles arose via relaxation processes in the complex 
systems of which the particles were once partfl 




4 An Estimate for the Relaxation Timescale 

It is important first of all to have a rough, order-of-magnitude estimate for the 
timescale over which relaxation to quantum equilibrium takes place. What we 
have in mind here is not the time taken to reach equilibrium but the timescale 
over which there is a significant approach to equilibrium.^ 

We may define a relaxation timescale r in terms of the rate of decrease of 
H(t) near t = 0. Because (dH/dt) Q vanishes one has to consider (d 2 H/dt 2 ) . 
Thus we define r by 1/t 2 = - (d 2 H/dt 2 ) /H (Valentini 1992), where (9) 
gives [d 2 H/dt 2 ) Q in terms of the initial state. Expanding Xq ■ V/o in a Taylor 
series within each coarse-graining cell of volume SV = (Sx) 3N , it is found that 



(d 2 H/dt 2 ) = -I(Sx) 2 /12 + 0{(Sx) i ) where/ - J dX (|* |7/o)|V(AV V/ )| 2 
(Valentini 2001a). Thus 





(10) 



For 8x small compared to the lengthscale over which Xq ■ V/o varies, r oc 1/Sx. 



5 For further details see Valentini (1991a,b; 1992; 2001a). 
6 Cf. the scattering time of classical kinetic theory. 
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Taking Ho <~ 1 (a mild disequilibrium) and crudely estimating I one obtains 

1 mh 2 



Sx (AP f 



(11) 



where AP is the (quantum) momentum spread of (Valentini 2001a). Given 
t oc 1/Sx this formula may also be obtained on dimensional grounds. Note that 
this result is merely a crude estimate. 

5 Numerical Simulations 

A numerical simulation of relaxation has been performed for the simplest pos- 
sible case of an ensemble of independent particles in a one-dimensional box 
(Valentini 1992). The results were surprisingly good, given that the particles 
cannot move past each other. This simple model is of course an unrealistic 
setting for relaxation: one should really consider an ensemble of complicated 
systems with many degrees of freedom. Nevertheless the model is instructive. 
Consider, then, a one-dimensional box with infinite barriers at x = 0, L 

and energy eigenfunctions </>„(x) — ^2/Lsm (mrx/L) (n = 1,2,3, ) with 

eigenvalues E n = ^(nn/L) 2 (units m = H = 1). The box contains an ensemble 
of independent particles, each guided by the same wavefunction ip(x, t). At t = 
the wavefunction was taken to be a superposition of the first M eigenfunctions, 
with amplitudes of equal modulus but randomly-chosen phases 9 n : 

M l 

i>o(x) = y^—=<j>n{x) exp(«6>„) 



, /M" 

n=l v 

The particles were taken to be uniformly distributed at t = 0: thus, po{x) = 1/L 
and of course po(x) ^ I^Ao^)! 2 - These initial conditions are sketched in Fig. 1, 
for the case M = 10 and L = 100. 

By numerical calculation of the trajectories (given by x — (d/dx) lining), 
the distribution p(x, t) at later times may be determined, while the Schrodinger 
evolution of ip(x, t) is just 



M 1 

i>{x,t) = V^=0„(x)expi(6>„ - E n t) 
■ 'M 



n=l 



For L — 100, ip(x,t) is periodic in time with period T P = 2n/Ei « 12,700 in 
our units. Because the particles cannot move past each other, each trajectory 
must recur with period Tp (to ensure that the equilibrium distribution recurs) 
and so any po{x) ^ \ipo{ x )\ 2 recurs as well. While this simple model exhibits 
a strong form of recurrence, nevertheless for times smaller than the recurrence 
time the approach to equilibrium is significant. 

An example is shown in Fig. 2, again for ten modes (M = 10). At t = 120, 
p has developed sharp peaks which coincide with the (smooth) maxima of \^\ 2 . 
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An experimenter with a 'blunt' measuring device, say with resolution Sx = 10, 
would conclude that p and |-0| 2 are roughly equal. The approach to equilibrium 
(at t <C Tp) is surprisingly good for such a simple, highly constrained modelj^] 

The evolution of H(t) has also been calculated (with C. Dewdney and B. 
Schwentker). An example is shown in Fig. 3, once again for M = 10, where H 
is defined with a coarse-graining length Sx — 1. The curve H = Hit) shows a 
strict decrease soon after t — (as it must); thereafter it decreases steadily but 
not monotonically, for t <C Tp. (Of course, as the recurrence time is approached, 
H necessarily increases as the initial conditions are restored.) 

It is clear from Fig. 3 that H initially decreases on a timescalc of order 
~ 100. This result compares well with our crude estimate (11). For here APq ps 
in/VS)M/L and with L = 100 and Sx = 1 our prediction is t ~ 2 x 10 5 /M 3 . 
For M = 10 we have r ~ 200 - in fair agreement with the numerical calculation. 

The timescalc over which H decreases - which may be quantified in terms 
of the time t$% taken for H to decrease by 5% - has been calculated (with C. 

Dewdney and B. Schwentker) for different numbers of modes, M = 10, 15, , 

40 (with fixed Sx — 1), and for different coarse-graining lengths ranging from 
Sx = 0.2 up to Sx — 2 (with fixed M = 20). The results agree quite well with 
our predictions that t 5 % cx 1/M 3 and t 5 % oc 1/Sx (Dewdney, Schwentker and 
Valentini 2001; Valentini 2001a).g 

The above model serves an illustrative purpose only. The particles cannot 
move past each other, and the approach to equilibrium is limited. Higher- 
dimensional models should show a much better approach to equilibrium. 

Calculations analogous to the above are currently being performed by Dewd- 
ney, Schwentker and Valentini for a two-dimensional box. Bearing in mind the 
generally chaotic nature of the trajectories for this system (Frisk 1997), an effi- 
cient relaxation to equilibrium is to be expected. 



6 Comments on 'Typicality' 

A completely different approach to explaining quantum equilibrium is based on 
a notion of 'typicality' for the initial configuration Xq mv of the whole universe 
(Diirr et al. 1992a, b). Here, if is the initial wavefunction of the universe, 

|^,unzD|2 - g t a ]j eil t be the 'natural measure' on the space of initial configura- 
tions. It is shown that quantum distributions for measurements on subsystems 
are obtained for 'almost all' Xq"™, with respect to the measure I^q™ 1 "! 2 - 

While Diirr et al. give a general proof, their approach may be illustrated by 
the case of a universe consisting of an ensemble of n independent subsystems 

7 For further details see Valentini (1992; 2001a). Note that here the quantum timescale for 
the evolution of \ip\ 2 is At = h/AE ~ 70, where AE is the quantum energy spread. 

8 In the first case, upon averaging over many trials with different random initial phases, a 
plot of ln<5% versus InM was found to have a slope of —2.6, compared with the predicted 
value of —3. In the second case, there are departures from linearity for Sx > 0.5 due to the 
large variation in the velocity field over some coarse-graining cells, the result r oc l/5x being 
valid only for small 5x. 
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(which could be complicated many-body systems, or perhaps just single parti- 
cles), each with wavefunction i)jo{x). Writing ^ niv = i/'o (^1)^0(^2) ipo(xn) 

and Xq uiv = (xi, x 2 , £3, , x n ), a choice of X^ nm determines - for large n 

- a distribution po(x) which may or may not equal |-0o(2O| 2 - 

Now it is true that, with respect to the measure |\J/g Tm ' | 2 ; as n — > 00 almost all 
configurations Xq"™ yield equilibrium po = \ipo\ 2 for the subsystems. It might 
then be argued that, as n — » 00, disequilibrium configurations occupy a vanish- 
ingly small volume of configuration space and are therefore intrinsically unlikely 
However, for the above case, with respect to the measure |\p™™| 4 almost all 
configurations X^ ntv correspond to the disequilibrium distribution po = |"0o| 4 - 
This has led to charges of circularity: that an equilibrium probability density 
|^r«mu|2 j s m e ff ec ^ being assumed for Xq"™; that the approach amounts to in- 
serting quantum noise into the initial conditions themselves (Valentini 1996) |] 




However, there is a certain rationale to Durr et aVs approach. It has been 
argued (S. Goldstein, private communication) that in statistical mechanics one 
always rules out unwanted initial conditions in some way, usually by deeming 
them exceptional (or 'untypical') with respect to a specified measure. This is 
perfectly true. But in the author's opinion, a class of initial conditions should 
be selected on empirical grounds. Even today we know only that, for all sub- 
systems probed so far, p = \ip\ 2 ± e where e is some experimental accuracy. 
There is certainly no reliable empirical evidence that p = \ip\ 2 exactly - or even 
approximately - in the very early universe, near the big bangj^ 

Note also that pilot-wave theory in exact equilibrium may never be suscep- 
tible to an experimental test of the basic law of motion (the de Broglie guidance 
equation); details of the trajectories may well be forever hidden from us, unless 
the human brain has an unexpected sensitivity to hidden variables (Valentini 
2001a). Restricting pilot-wave theory to quantum equilibrium is as artificial as 
it would be to restrict classical mechanics to thermal equilibrium; and a huge 
amount of potentially new physics is lost thereby. 

Finally, Durr et al. focus on initial configurations for the whole universe 
while in Sect. 2 we emphasised initial distributions for subsystems. One might 
think that the former description is more correct because in the early universe 
there are no independent subsystems (that is, subsystems with their own wave- 
function). But even leaving aside the point that the asymptotic freedom of 
interactions at high energies allows one to treat particles at very early times as 
essentially free, the two modes of description are really interchangeable: even 
if there are no independent subsystems, one can imagine (theoretically) what 

9 Note that if the word 'typicality' is replaced by 'probability', the result of Durr et al. 
(1992a, b) becomes equivalent to the 'nesting' property proved by Valentini (1991a), which 
states that an equilibrium probability for a many-body system implies equilibrium probabili- 
ties for extracted subsystems, as briefly discussed above. 

10 In fact, the horizon problem in cosmology (the uniformity of the microwave background 
over regions larger than the classical causal horizons of a Friedmann universe) suggests that 
noncquilibrium p may have existed at early times, the resulting nonlocality playing a 

role in homogenising the early universe (Valentini 1991b; 1992; 2001a). Inflation attempts to 
remove the horizon problem by changing the early expansion; but as yet there is no satisfactory 
inflationary theory. Early quantum disequilibrium offers an alternative approach. 
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would happen if (say) single particles were extracted from a gas of interact- 
ing particles, prepared with a certain wavefunction ip, and had their positions 
measured. Whether the measured distribution p is equal to \tp\ 2 or not will 
depend on the configuration of the whole system prior to the extraction. By 
this (purely theoretical) device, one could characterise the total configuration as 
'equilibrium' or 'nonequilibrium' (Valentini 2001a). Our Sect. 2 could then be 
rephrased in terms of configurations rather than distributions: our 'empirical' 
approach would amount to trying to find out what sort of initial configuration 
could explain the statistics observed today - and again, while an initial 'equi- 
librium' configuration is a possibility, there are clearly many 'nonequilibrium' 
configurations at t — that explain the observations just as well. 



7 The Early Universe. Suppression of Relax- 
ation at Early Times 

At first sight one might think that, even if the universe did begin in quantum 
nonequilibrium, relaxation p — > \%l)\ 2 would take place so quickly (in the extreme 
violence of early times) as to erase all record of nonequilibrium ever having 
existed. However, closer analysis reveals a plausible scenario whereby quantum 
nonequilibrium could survive to the present day, for some types of particle left 
over from the big bang, leading to observable violations of quantum theory. 

The scenario has three stages: (i) depending on how the relaxation timescale 
r scales with temperature T, it is possible that relaxation will be suppressed at 
very early times by the very rapid expansion of space, in which case (ii) those 
particles that decouple at very early times will still be out of quantum equilib- 
rium at the time of decoupling, after which (iii) the free spreading of particle 
wavefunctions stretches any residual nonequilibrium up to larger lengthscales 
that are within reach of experimental test. This scenario will now be briefly 
sketched, in this section and the next]"] 

In the standard hot big bang model, the early universe contains what is essen- 
tially an ideal gas of effectively massless (relativistic) particles. It is instructive 
to recall the standard reasoning about how thermal equilibrium between distinct 
particle species is achieved: the interaction timescale U nt - that is, the mean 
free time between collisions - must be smaller than the expansion timescale 
t exp = a/a ~ (1 sec)(l MeV/kT) 2 , where a(t) cx t 1 ! 2 is the scale factor and 
t is the time since the beginning.^ As is well known, this condition need not 
always be satisfied: for instance, if the cross section a(T) falls off faster than 
1/T at high temperatures then ti nt ~ 1/nca - with a number density n ~ T 3 
(in natural units where H = c = 1) - will fall off slower than 1/T 2 so that, at suf- 
ficiently high temperatures, U n t > t exp and the particles do not thermalise. The 
functions a — o~(T) are determined by particle physics; once they are known, it 
is possible to determine the thermal history of the universe, and in particular 

11 For further details see Valentini (2001a, c). 
12 See any standard cosmology text. 
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to deduce when particle species fall out of thermal contact with each other. 

Similar reasoning applies to the relaxation p — > IV'I 2 - There will be a 
temperature-dependent timescale r = r(T) over which relaxation takes place. 
There will also be a competing effect due to the expansion of space: as space 
expands wavefunctions are stretched, and it is easy to see in pilot-wave theory 
that this results in a proportionate expansion of the disequilibrium lengthscale. 

To see this, consider an initial distribution po(x) with disequilibrium on a 
small lengthscale S - We may write po{x) = \ipo(x)\ 2 fo(x), where ipo{x) has 
some width A and fo(x) ^ 1 varies over distances So <C A (so that coarse- 
graining over distances much larger than Sq yields an equilibrium distribution) . 
If the wavefunction expands up to a width A(t), then because / is conserved 
along trajectories, and because particles initially separated by a distance fen 
are later separated roughly by 5x(t) ~ (A(t)/A ) Sxq, it follows that deviations 
/ 1 occur on an expanded lengthscale 5(t) ~ (A(t)/An)5n- This is true 
irrespective of whether the wavefunction expands because of the expansion of 
space or simply because of its own natural free spreading over time.0 

Thus relaxation p — > \tp\ occurs if r < t exp , while it is suppressed if r > t exp . 
The issue now depends on how r scales with T. Unfortunately, there are as yet 
no reliable calculations that can tell us this. But we can make a crude estimate. 

Given t oc 1/Sx, for massless particles r may be estimated on purely di- 
mensional grounds to be t ~ (l/Sx)H 2 /c(Ap) 2 where Ap is the particle mo- 
mentum spread. Taking Ap ~ kT/c we have r ~ (\/5x)h 2 c/{kT) 2 which is 
larger than t exp ~ (1 sec)(l MeV/fcT) 2 for all T if Sx < 10l P (where l P is 
the Planck length). However, it would be more realistic to apply our estimate 
for t to a coarse-graining length 5x that is relevant to the energetic processes 
at temperature T, the natural choice being the thermal de Broglie wavelength 
8x ~ Hc/kT (the typical width of particle wavepackets). We then have r ~ H/kT 
- a plausible result, according to which relaxation is suppressed (r > t exp ) if 
kT > 10 18 GeV« O.lkTp or t < 10t P (where T P and t P are the Planck tem- 
perature and time). Our simple estimate suggests that relaxation to quantum 
equilibrium began one order of magnitude below the Planck temperature. 

Note that in the standard Friedman expansion considered here, the rate 
at which space is expanding becomes infinite as t — > 0. Of course, our 
relaxation rate 1/t oc T also tends to infinity as t ~ * (and T ~ * oo), but this 
is offset by the expansion rate l/t exp oc T 2 which tends to infinity even faster. 
Note also that while our estimate is very crude and should not be taken too 
seriously, it does illustrate the key point that, if r decreases with temperature 
more slowly than 1/T 2 , then relaxation will be suppressed at very early times. 

It is therefore plausible that any particles that decouple soon after the Planck 
era will not have had time to reach quantum equilibrium, the extreme violence 
of that era being offset by the even more extreme expansion of space. And at 
the time of decoupling, such particles are expected to show deviations p ^ \ip\ 2 
on a lengthscale 5x ~ hc/kT P = l P 10~ 33 cm. 

13 This is consistent with the //-theorem, for H measures negentropy relative to \ip\ 2 : if \ip\ 2 
itself expands, the absolute disequilibrium lengthscale can increase without increasing H. 
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8 Residual Disequilibrium Today. Experimen- 
tal Tests 



After decoupling, particle wavefunctions undergo a huge expansion - partly 
due to the expansion of space itself and partly due to the free spreading of 
the wavepacket through space. This results in a proportionate stretching of 
disequilibrium to much larger lengthscales, as shown above.0 

For example, relic gravitons are believed to decouple at kTd ec ~ 10 19 GeV 
or at redshift Zd ec ~ Td e c/T now ~ 10 32 . The subsequent expansion of space 
alone, by a linear factor of ~ 10 32 , will stretch the disequilibrium lengthscale 
up to ~ 1 mm. Of course, there seems to be little hope of detecting the ~ 
1° K graviton background directly in the near future - still less of testing it 
for violations of quantum theory. However, it is expected that there are other, 
more exotic particles that decoupled soon after Tp. (Supersymmetry and string 
theory predict a plethora of new particles at high energies.) These may be no 
easier to detect directly. However, some of them might decay at later times 
into more easily detectable particles such as photons, perhaps by annihilation 
X + X — > 27 or by a supersymmetric decay X — > 7 + 7 into a photon and a 
photino. Nonequilibrium for the parent particles should yield nonequilibrium 
for the decay products as well. Thus, we would suggest testing the Born rule p = 
\tj}\ 2 for photons produced by the decay of exotic relic particles from the Planck 
era. For example, in a two-slit interference experiment, such photons might 
produce an interference pattern that deviates from the quantum prediction. 

Experiments are under way searching for exotic relic particles supposed to 
make up the 'dark matter' pervading the universe. Some of these experiments 
involve searching for decay photons. However, the usual dark matter candidates 
(such as neutrinos, neutralinos, axions or gravitinos) are expected to decouple 
much later than tp. Nevertheless, our knowledge of particle physics beyond 
the standard model is so uncertain that, for all we know, there might exist 
an appropriate relic particle X - that decouples soon after tp and partially 
accounts for dark matter. If such particles were detected, they or their decay 
products would be our prime candidates for particles violating the Born rule. 

In the meantime one might also consider the relic photons that make up the 
microwave background. But why consider particles that decoupled at t ~ 10 yr 

- much later than tp - when our estimates suggest that any disequilibrium will 
have been erased? One answer is that the hugely expanded wavepackets of 
relic particles provide an interesting test of the Born rule in extreme conditions. 
On general grounds (irrespective of the suggestions made here about hidden 
variables and early nonequilibrium), it would be worth testing quantum theory 
in such unusual circumstances, where quantum probabilities have spread over 

14 Decoupling is of course never exact. Because of small residual interactions at all times, 
the wavefunctions of relic particles contain tiny scattering terms that perturb the trajectories 

- possibly 're-mixing' the expanding disequilibrium. Calculation of the effect of scattering 
by a tenuous medium shows that re-mixing does not happen: the trajectory perturbations 
grow only as t 1 / 2 and cannot overcome the linear (oc t) free expansion of the disequilibrium 
lengthscale (Valentini 2001a, c). We may therefore safely ignore such residual interactions. 
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intergalactic distances. Relic photons decouple at kTd ec ~ 1 eV (zdec ~ 10 3 ) and 
at decoupling their wavepackets have widths of order ~ fic/(l eV)^ cm. 
If the packets have (as is widely assumed) spread essentially freely at the speed 
of light for ~ 10 10 years then today they will have a width ~ 10 28 cm and so 
the packets will have expanded by ~ 10 33 (ignoring the relatively small effect 
of the expansion of space by a factor ~ 10 3 ). It can be argued - on general 
grounds, independent of pilot-wave theory - that for cosmological microwave 
photons the quantum probability today on lengthscales ~ 1 cm could contain 
traces of corrections to the Born rule which may have existed (for whatever 
reason) on the Planck scale at the time of decoupling (Valentini 2001a, c).^ 

As suggested elsewhere (Valentini 1996), it would then be worthwhile to test 
quantum theory for photons from the microwave background. In a two-slit in- 
terference experiment with single relic photons, the usual quantum interference 
pattern - given by p = \i/j\ 2 — might show an anomalous blurring^ 

Generally, corrections to the Born rule for any photons from deep space 
would produce a number of observable anomalies. In particular, the functioning 
of some astronomical instruments might be affected: diffraction-grating spec- 
trometers could produce unreliable readings, and the diffraction-limited perfor- 
mance of some telescopes might be impaired]^ 



9 Outlook 

It is clear that much remains to be done to develop the above ideas fully. Other 
approaches to subquantum statistical mechanics (based for example on external 
perturbations) remain to be developed. Properties of the trajectories such as 
ergodicity and mixing (in a rigorous sense) should be investigated. 

The process of relaxation to quantum equilibrium should also be studied 
further, in particular in the early universe. It remains to be seen if the relaxation 
timescale r really does decrease with temperature more slowly than 1/T 2 , so 
that relaxation is suppressed at very early times. An experimental test of the 
Born rule p = \ip\ 2 for relic cosmological particles seems feasible, in particular 
for photons from the microwave background (Valentini 2001a,c). 

As for the issue of chance in physics, the central conclusion of this work is 
that, like the cosmic microwave background, quantum noise is a remnant of the 
big bang. And just as the microwave background has been found to have small 
nonuniformities in temperature, so the 'quantum background' - the quantum 
noise that pervades our universe - may have small deviations from the Born 
rule p = \ip\ 2 , and should be probed experimentally. 

15 Our point here is that the line of argument given in this paper has led us to propose 
experimental tests that are actually worthwhile in their own right, because they would probe 
quantum theory in new and extreme conditions. 

16 Incoming photons will have very nearly plane wavepackets, so there is no ambiguity as 
to what quantum theory would predict in such an experiment. Note also that it has recently 
become possible to detect single photons in the far-infrared region (Komiyama et al. 2000). 

17 For details see Valentini (2001a). 
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FIGURE CAPTIONS 

Fig. 1. Initial conditions at t = 0. We have plotted \ip\ 2 for a superposition 
of ten modes with amplitudes of equal modulus and random phase. The initial 
density p is taken to be uniform. 

Fig. 2. At t = 120. There is a strong coincidence between the sharp 
peaks of p and the smooth peaks of \il>\ 2 , indicating an approximate approach 
to equilibrium on a coarse-grained level. 

Fig. 3. Plot of the coarse-grained ff-function against time. 
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